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Formal solutions to the equations of motion of the one-dimensional, unsteady flow of a viscous, 
compressible, heat-conducting gas are presented. The relationship between the existence of these 
solutions and the problem of hydrodynamic stability is discussed. 


INTRODUCTION 

T HE present investigation aims to demonstrate 
that the one-dimensional flow of a viscous, 
compressible, and heat-conducting fluid can be 
represented formally by the heat-conduction equa- 
tion in the stream function-time plane. No explicit 
solutions are obtained; nevertheless the behavior 
of the fluid can be deduced without actually solving 
the equations of motion. 

In order to simplify the analysis, the stream 
function is introduced as a streamwise coordinate 
and the Prandtl number, Pr, is taken to be f. This 
coordinate transformation is patterned after the 
von Mises transformation. 1 The momentum equa- 
tion is nonlinear in the space-time plane by virtue 
of its convective acceleration term; in the stream 
function-time plane the term drops out and the 
equation assumes the form of the well-known 
diffusion equation in a homogeneous medium. This 
equation is, of course, linear, and the classical 
methods of dealing with this type of equation can 
be brought to bear on the investigation. 


ANALYSIS ' 


The equation expressing the conservation of 
energy 2 is given by 


fj w + w 
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( 1 ) 


where c v is the specific heat at constant pressure, X 
is the thermal conductivity, and v is the streamwise 
component of velocity. Substituting Pr \/c p for p, 
and setting Pr = f yields 



1 R. von Mises, Mathematical Theory of Compressible 
Fluid Flow , (Academic Press, Inc., New York, 1958), pp. 157- 
158. 

1 J. 0. Hirschfelder, and C. F. Curtiss, J. Chem. Phys. 17, 
1076 (1949). 


where H = c p T + | v 2 . The momentum equation 
for Pr = f is simply 
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The continuity equation is given by 
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The stream function, defined by 

p = d^/dx, pv = -df/dt, (5) 

may be introduced as a streamwise coordinate, 
yielding 
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Application of this transformation to the energy 
equation gives 
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( 6 ) 


while the momentum equation reads 
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• where k 2 = \p/c v . 

It is to be noted that the momentum equation 
(which is coupled to the energy equation only 
through the pressure gradient term) is now linear 
while the energy equation remains nonlinear. Our 
plan is to treat the pressure gradient term in the 
momentum equation as a known function of \f/ and 
t , and to treat Eq. (7) as if it were a nonhomogeneous 
heat conduction equation. We shall be able to 
demonstrate that without actually knowing the 
dependence of p on \f/ and t , and by imposing only 
certain physically plausible restrictions, we can 
establish whether or not the equations of motion 
possess solutions that are bounded for all times. 

The linear nature of Eq. (7), together with the 
hypothesis concerning the pressure gradient, per- 
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mits the use of the principle of superposition. Thus 
we may first solve the homogeneous equation with 
inhomogeneous boundary conditions and then solve 
the inhomogeneous equation with homogeneous 
boundary conditions. The complete solution is a 
linear combination of the two. Calling the solution 
to the homogeneous equation v u we obtain by the 
method of Green's function 
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dr, ( 8 ) 


where v x = at t = 0. We are, of course, 

interested in the behavior of the solution as t *— > co ; 
thus we have to take the limit of the integral. Since 
the integration is independent of t , and, moreover, 
since the integrand is uniformly continuous with 
respect to t , it is permissible to interchange the order 
of integrating and of taking the limit, provided the 
improper integral (8) converges. A sufficient con- 
dition for the convergence of the integral is the 
convergence of 


V 2 = exp ( — k 2 s 2 t) J exp (k 2 s 2 r)sP(s, r) dr. (12) 

Taking cognizance of the fact that we are dealing 
Math a real as opposed to an ideal gas, we now 
assert on a physical basis that whatever the pressure 
may be — and for that matter, its image function 
P(s, t ) — it must be bounded. The rationale for this 
is simply that a viscous and heat conducting gas 
cannot support an infinite pressure. Setting the 
maximum value of the pressure equal to M 2 permits 
us to perform the integration and obtain 


M 2 s[l — exp ( — k 2 s 2 t)] 

kV 


(13) 


Before performing the inversion integration to 
return to the physical plane, we let t — » °° t and thus 
obtain 


lim v 2 < (M//c) 2 . 


(14) 


f im dr. (9) 

In other words, any initial velocity distribution, 
compatible with the boundary conditions, that 
satisfies the convergence condition is a permissible 
one. In turn, any permissible initial distribution 
satisfying the . boundary conditions and the con- 
vergence criterion will die out as thus 

lim v x — 0. ■ (10) 

- — ■ • * — oo . 

The particular solution to the inhomogeneous 
equation v 2 (\j/ 1 t) is obtained with the aid of the 
Laplace transform. We take the Laplace transform 
of Eq. (7) with respect to \p t calling the transform 
variable s and the transformed functions by their 
respective capital letters, V 2 and P. We thus obtain 
for the transformed momentum equation 


This inequality asserts that v 2 is bounded. Our 
apparent success in dealing with the momentum 
equation by showing the existence of a bounded 
solution obviates the need to consider the energy 
equation, which can be dealt with in precisely the 
same fashion. 

DISCUSSION 

We have demonstrated that the unsteady equa- 
tions of motion possess only bounded solutions. A 
bounded solution is also a stable one if after the 
initial conditions are changed slightly the new 
solution will approach the original one asymp- 
totically. Since the initial conditions are involved 
only in the solution to the homogeneous equation, 
which we have shown to be zero * provided that 
f(r)dr converges, it follows that the solution is 
stable. 


— + kW 2 = sP(s, t), (11) 
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